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Excitation of a metallic ferromagnet such as Ni with an intensive femtosecond laser pulse causes 
an ultrafast demagnetization within approximately 300 fs. It was proposed that the ultrafast de¬ 
magnetization measured in femtosecond magneto-optical experiments could be due to relativistic 
light-induced processes. We perform an ab initio investigation of the influence of relativistic effects 
on the magneto-optical response of Ni. To this end, we develop, first, a response theory formula¬ 
tion of the additional appearing ultra-relativistic terms in the Foldy-Wouthuysen transformed Dirac 
Hamiltonian due to the electromagnetic field, and, second, compute the influence of relativistic 
light-induced spin-flip transitions on the magneto-optics. Our ab initio calculations of relativistic 
spin-flip optical excitations predict that these can give only a very small contribution (< 0.1%) to 
the laser-induced magnetization change in Ni. 

PACS numbers: 71.15.Rf,78.47.J-,78.20.Ls 


I. INTRODUCTION 

Ultrafast laser-induced demagnetization of metallic 
ferromagnets was discovered in 1996 by Beaurepaire et 
al. [1], who observed that a ferromagnetic Ni film could 
be demagnetized to ^^50% in about 300 femtoseconds af¬ 
ter excitation with a short laser pulse. This surprising 
discovery was followed by many pump-probe magneto¬ 
optical experiments on elemental metallic ferromagnets 
that confirmed the phenomenon of laser-induced demag¬ 
netization (see e.g., Refs. [2-8]). More recently, ultrafast 
laser-induced demagnetization has been studied in multi¬ 
sublattice or multilayer materials employing element re¬ 
solved probing techniques in the extreme ultraviolet and 
soft x-ray regimes [9-14]. 

The discovery of ultrafast laser-induced demagnetiza¬ 
tion led to an intensive debate on what the underly¬ 
ing microscopic mechanism of the ultrafast dissipation 
of spin angular momentum could be [15-17]. Several 
mechanisms have been proposed to explain the ultra¬ 
fast demagnetization and these continue to be discussed 
[18-22]. An early microscopic explanation was based 
on direct transfer of angular momentum from the light 
involving the spin-orbit interaction [18]. Another pro¬ 
posed mechanism is spin dissipation through fast Elliott- 
Yafet electron-phonon spin-flip scatterings [20]. Other 
proposals are electron-magnon spin-flip scattering [19] or 
electron-electron spin-flip scattering [21]. A different sce¬ 
nario is based on the laser-generation of superdiffusive 
spin currents that transport spin angular momentum out 
of the excited ferromagnetic film, thus reducing its net 
magnetization [22, 23]. Other explanations have focused 
on the direct action of the laser on the electron’s spin, 
causing either a direct, laser-induced spin-flip [24] or a 
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change of the spin through an ultra-relativistic spin-light 
interaction [25]. 

Despite the still ongoing debate on the mechanism of 
ultrafast demagnetization it has been shown that spin dy¬ 
namics simulations within the Landau-Lifshitz-Gilbert, 
Landau-Lifshitz-Bloch, or Landau-Lifshitz-Baryakhtar 
formulations [26-29] can be used to describe ultrafast 
laser-induced demagnetization in alloys when a suffi¬ 
ciently large and fast dissipation of spin angular momen¬ 
tum is assumed. 

To establish accurately how much demagnetization can 
be caused by one of the aforementioned mechanisms 
density-functional theory (DFT) based electronic struc¬ 
ture calculations are indispensable. Recently, DFT-based 
investigations have been performed for the Elliott-Yafet 
electron-phonon spin-flip scattering in transition metal 
ferromagnets [30-33]. The ab initio calculations pre¬ 
dicted relatively small demagnetization rates; this gave 
rise to modified proposals, in which in addition an ul¬ 
trafast reduction of the exchange splitting needed to be 
taken into account to explain the observed demagnetiza¬ 
tion [34, 35]. Another recent computational investigation 
suggested that a combination of spin-flip electron-phonon 
and electron-magnon scatterings could explain the mea¬ 
sured demagnetizations [36]. 

A demagnetization scenario involving the direct, rela¬ 
tivistic spin-photon interaction was proposed a few years 
ago [25]. In this proposal ultra-relativistic terms stem¬ 
ming from the Dirac Hamiltonian provide a coupling be¬ 
tween the electromagnetic field of the pump laser pulse 
and the spins of electrons in the material [25, 37, 38]. 
Model calculations of this mechanism have recently been 
made for transitions from the 2 s to 2p levels of a hydrogen 
atom [39]. A full ah initio investigation of the influence 
of the relativistic spin-photon interaction on the magne¬ 
tization and magneto-optical response has not yet been 
made. 

Here, we report an ab initio investigation of the in- 
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fluence of the relativistic spin-photon interaction. We 
present first analytic theory to analyze which terms are 
the relativistic terms that are involved in the coupling 
of the spin and photon fields. A notable difference as 
compared to other recent investigations [25, 39] is the di¬ 
rect consideration of the exchange field in our approach. 
Also, we investigate the influence of the additional rel¬ 
ativistic terms on the response theory equations for the 
magneto-optical spectrum. The derived expressions are 
employed in ab initio calculations of the magneto-optical 
Kerr effect (MOKE) of Ni. Our calculations underline 
that the influence of relativistic, laser-induced spin-flips 
is present, but is quite small and can thus not account 
for the substantial amount demagnetization that is ob¬ 
served in femtosecond pump-probe magneto-optical mea¬ 
surements. 

In the following we first provide a derivation of the rel¬ 
ativistic spin-photon interaction starting from the Dirac 
equation (Sec. IIA and IIB). In Sec. IIC the correspond¬ 
ing momentum operator is derived and results of ab initio 
calculations for Ni are presented in Sec. Ill, and conclu¬ 
sions in Sec. IV. 


II. THEORY 


A. The Dirac-Kohn-Sham equation 


To include the relativistic light-spin coupling effects in 
the calculation of the magneto-optical Kerr spectra we 
consider the Dirac-Kohn-Sham (DKS) equation [40-42] 

H\^) = {E-mc^m ( 1 ) 

with H being the DKS Hamiltonian, 

i/= ca ■ p-f (/3 — 1) TOc^-I-V-|-/is/SS • B™. (2) 


Here S = 1 ®(t is the spin operator in the Dirac bi-spinor 
space, V is the unpolarized Kohn-Sham selfconsistent po¬ 
tential, is the spin-polarized part of the exchange- 
correlation potential in the material, p = —ih^, 1 is 
the 4x4 identity matrix, and pB is the Bohr magneton, 
Pb = The matrices 


a = 





0 
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momentum. It is thus not a proper magnetic field and 
cannot be represented by a vector potential [41[. There¬ 
fore it is not included as a vector potential in the lin¬ 
ear momentum, i.e. p — eA^^. However, to account for an 
external electromagnetic perturbation (for instance being 
due to a laser source) we introduce the vector potential 
A(r, t), leading to 

H = coL-{p-eA)^{(3- 1) mc^ + V + pb/3S• 5=^=. (3) 

In Appendix A we show that, in the nonrelativistic limit, 
this form of the DKS equation leads to a Hamiltonian 
where the external magnetic field B{r,t) = V x A{r,t) 
couples to the spin S (= |(t) and orbital angular mo¬ 
mentum L operators, but the exchange field B^"^ couples 
only to the spin operator. 

The aim of this work is to investigate the influence of 
relativistic terms - that lead to a spin-photon field cou¬ 
pling - on the MOKE spectra. To elucidate the terms 
that involve both spin degrees of freedom and the ex¬ 
ternal electromagnetic field we rewrite the Hamiltonian 
in Eq. (3), as a semirelativistic expansion in terms of 
order of 1/c^. Such rewriting of the DKS equation can 
be achieved in two ways. The small component of the 
wavefunction \tjj-) can be eliminated exactly, leading to 
an equation which is fully equivalent to the DKS equa¬ 
tion, but for the large component ['0+) only [43[. This 
equation can subsequently be expanded in orders of 1 jt ?. 
Alternatively, one can use the Foldy-Wouthuysen (FW) 
transformation approach [42, 44[, however, it needs to 
be extended to the case where an exchange field Bi^'^ is 
present, which was not done before. We note that apart 
from the exact transformation of Kraft et al. [43[, a Green 
function technique was applied by Crepieux and Bruno 
[45[ to obtain the semirelativistic Hamiltonian. However, 
they did not start from the DKS Hamiltonian as given 
in Eq. (3), but instead added the external magnetic field 
to the exchange field and did not have a vector potential 
in the momentum, p — eA (but only p). As discussed 
further below, they obtained several similar terms, yet 
not all that follow from the FW transformation. 

In the following we employ the Foldy-Wouthuysen 
transformation to derive the Hamiltonian terms that give 
rise to a spin-photon field coupling. 


B. The time-dependent Foldy-Wouthuysen 
transformation 


are the well known Dirac matrices, with cr the Pauli spin 
matrices and 1 the 2x2 identity matrix. The fully rela¬ 
tivistic state \ijj) in Eq. (1) is the Dirac bi-spinor. 



An important point to observe is that in the DKS equa¬ 
tion the exchange field B^^ is different from the standard 
magnetic field, as it obviously acts only on the spin de¬ 
gree of freedom and does not couple to the orbital angular 


To make a clear distinction between pure non¬ 
relativistic Schrodinger-like terms and relativistic terms 
(up to the order of 1/c^) we use the Foldy-Wouthuysen 
transformation [42[. Formally, the time-dependent FW 
transformation can be expressed as 

Hpw = (^H - ih^^ (4) 

where is a unitary operator that transforms the 

DKS Hamiltonian to a block diagonal form, where each 
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block is 2 X 2. The four-component Dirac Hamiltonian 
is diagonalized under the assumption that, at all points 
in the conhguration space, two of the spin components 
are much smaller than the other two. This assumption is 
valid if the kinetic and potential energy of the electron is 
much smaller than the rest mass energy of the electron. 
Since we are interested only in the “positive energy” so¬ 
lutions we will retain only the upper 2x2 component of 
the Hamiltonian (the large component of the Dirac bi¬ 
spinor). To find the transformed Hamiltonian in a 1/c^ 
expansion we write the DKS Hamiltonian as 

H = {(3-l)mc^ + 0 + S (5) 

with O = ca ■ {p — eA) as an odd operator (i.e., off- 
diagonal in the particle-antiparticle Hilbert space) and 
£ = V -|-p_b/ 35] - as an even operator (diagonal in 
the same space). Next, we consider the FW operator, 

Ufw = ~ 

To obtain an expansion in orders of 1 ji? we use a Taylor 
expansion of the operator ~ 1 ± HIfw + 0{l/c^). 


This leads to a transformed Hamiltonian 

H' = {(3- 1) mc^ + 0' + £', (7) 

where O' is the transformed odd part which is of the or¬ 
der of 1/c? and £' is the transformed even part. Repeat¬ 
ing this procedure two times with the operators = 
— 2 mc^ ^ further transformed H", O", and £" 

and then working with the operator f30" 

on the transformed Hamiltonian H” we get rid of the 
odd terms up to the order of 1/c®. After these transfor¬ 
mations we obtain a transformed Hamiltonian written in 
terms of the original odd and even parts 

+ ^-^APAO,S]+^^]. ( 8 ) 

Substituting the explicit form of the operators O and £, 
retaining only the terms up to the order 1/c^ and keeping 
in mind that, for the vector potential of the external elec¬ 
tromagnetic field, B = W X A and E = — we arrive 
at the Hamiltonian restricted to the large component of 
the Dirac bi-spinor, 


Hpw = 


{p - eAY 


2m 

i 

Mb 

inn?c^ 

ipB 

4m^c^ 


V - pb(t ■ B^''- ^ib(t ■ B- 
eh 


{p - eA) 


1 


eh^ 

(p^v) - V • E 


(T ■ [pV) X (p — eA) — 




Sm^c^ ^ 

cr ■ {E X [p — eA) — (p — eA) x Ej 


I [p^((T • -I- 2cr ■ {pB^Y ■ (p — eA) + 2{p ■ B^Y iP~ eA) -I- • [p — eA)] cr ■ [p — eA) | 

[(pxB-).(p-eA)]. (9) 


Note that, when the momentum operator and a (vector) 
function are enclosed in round brackets the momentum 
operator acts only on this function. This Hamiltonian is 
an extension to the conventional Pauli Hamiltonian (see 
Appendix A) yet, including all the Xjc? terms and all the 
terms involving B'^'^ to the same order. We note in ad¬ 
dition that the Hamiltonian (9) is quite different from 
the Hamiltonian given by Bigot et al. [25], as they did 
not consider the magnetic exchange interaction, which 
however is the strongest magnetic interaction in a ferro¬ 
magnetic material as Fe, Co, or Ni. In a further work 
Vonesch and Bigot [39] considered a static homogeneous 
applied magnetic field, expressed by a vector potential, 
as well as a time-varying vector potential to describe the 
electromagnetic field. Such static homogeneous magnetic 
field is nonetheless different from the exchange field B ^'', 
as the latter, as mentioned in the previous section, can¬ 
not be included by means of a vector potential [40, 41[. 
Specifically, compared to the extended Pauli Hamilto¬ 
nian (9), Vonesch and Bigot obtain the first, second, 


fourth, sixth, and eighth terms, as well as a term sim¬ 
ilar (but not identical) to the ninth term. As they didn’t 
include an exchange field, they did not obtain any of the 
terms containing B^‘^, but had a contribution cr • Bext 
due to the static applied homogeneous field. In addition 
they found an additional term, Ab ■ Aext, the product 
of the two vector potentials of the electromagnetic radia¬ 
tion and the constant external magnetic held. This term, 
which stems from writing out {p — e[AB + Aext])^, does 
not appear in our formulation where the exchange held 
is not represented by a vector potential. 

Our extended Pauli Hamiltonian (9) can further be 
compared with the Hamiltonian obtained by Crepieux 
and Bruno [45] using a Green’s function-based method. 
They considered two different Hamiltonians, one with an 
effective (including exchange) vector potential Aos, and 
the other one with an effective magnetic held, Be®. In 
the hrst case they obtained the same terms as we do from 
the introduction of the external vector potential, how¬ 
ever, the feasibility of the DFT-based formulation of such 
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Hamiltonian appears to be an open question. In the sec¬ 
ond case, they apply their method to an effective Hamil¬ 
tonian where the orbital magnetic effects are neglected, 
and make the critical assumption = V x Ag®; how¬ 
ever, on account of its nature, the exchange part of their 
BeB does not satisfy the Maxwell equations as a nor¬ 
mal magnetic field does (see e.g. Ref. [41]). In our case 
we used a vector potential to account for an external 
field without any need to use the mentioned critical as¬ 
sumption, thus giving in our Hamiltonian the proper in¬ 
dependent terms and terms linear in both, the external 
vector potential and B^^ fields, which are missing in the 
approach of Ref. [45]. 

The Hamiltonian (9) looks cumbersome at first sight 
but its physical content is readily explained. 

• The first and second terms comprise the usual 
Schrodinger Hamiltonian for a particle in an exter¬ 
nal field V (which in this case is a selfconsistent po¬ 
tential) and where the minimal coupling with an ex¬ 
ternal vector potential A{r,t) is present (this may 
represent, in general, any kind of external electro¬ 
magnetic field). 

• The third term is a Zeeman-like term due to the 
presence of the magnetic exchange field. 

• The fourth term is the standard Zeeman term with 
the external magnetic field. 

• The fifth term is the relativistic mass correction. 

• The sixth and the seventh terms are respectively 
the Darwin terms related with the selfconsistent 
potential V and the standard Darwin term arising 
from the external perturbation. 

• The eighth and the ninth terms are those which 
in a central potential V give rise to the spin-orbit 
coupling. 

• All the remaining terms, except the last, can be 
seen as corrections to the spin-orbit coupling due 
to the spin-polarized exchange field. This is more 
apparent using the identity 

V 

rp = rpr - X L, (10) 

r 

where pr = -ih-^ is the spatial part of the mo¬ 
mentum operator. 


• The last term depends on the B^"^ field but is in¬ 
dependent of the spin. 

The terms which involve a direct coupling of the spin 
to the external electromagnetic field are the fourth, 
eighth, ninth, and tenth ones. Via these terms it would 
in principle be possible to control the spin of electrons 
in a magnetic material by applying an external laser field. 

C. Strategy for the numerical implementation 

Our aim is to implement the above-derived relativis¬ 
tic terms in a suitable formalism for ab initio calcula¬ 
tions. An adequate way to achieve this is to consider the 
change of the Hamiltonian due to the applied electromag¬ 
netic field, which then is treated as a perturbation within 
Kubo linear-response theory to obtain the corresponding 
optical conductivity tensor (see, e.g. [46]). 

The standard strategy for the derivation of the con¬ 
ductivity tensor consists in gathering the external mag¬ 
netic vector potential related linear terms in an interac¬ 
tion Hamiltonian [47], 

S{Hi) = - Jdr j ■ dA, (11) 

and then to rewrite the current density operator j in 
terms of the momentum operator. This procedure is 
straightforward in the fully relativistic approach [see Eq. 
(3)] because the momentum operator, which is the con¬ 
jugate of the position operator, is given as 

2771 

IId = — —[r,H]=mca, (12) 

and the variation of the DKS Hamiltonian [Eq. (3)] re¬ 
sults easily 

SHi e , . 

J = —j^ = ecQ; = —IId. (13) 

oA m 

However, in the semirelativistic limit this equivalence 
breaks down. The expression of the conjugate momen¬ 
tum operator, obtained using the position-momentum 
conjugation relation starting from the extended Pauli 
Hamiltonian [Eq. (9)], is: 


lip =p-h 


1 

4mc^ 


+ i(TX {pV) + • (pB"") + (p • B^^)a- + 2B"^(cr • p) + 2 (t(B"" • p) + i{p x B^^)} 

ml J 

(14) 


where all the terms in the square brackets are due to rel- the square brackets are related to the relativistic mass 
ativistic corrections. The first and the second terms in correction and the spin-orbit coupling, respectively, and 
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can be obtained by means of the standard FW trans¬ 
formation in the absence of the exchange field; all the 
remaining terms are new and stem from relativistic cor¬ 
rections and the exchange field in the DKS equation. To 
reformulate this momentum operator, we use that it has 
been seen in Ref. [48] that the current density operator 
in the semirelativistic limit can be written in the form 

3{r) = j-K+jsp = —Hp - xer + crxp) (15) 

m 2m 

where ^ means that the operator p operates to the left. 
Here the second and third terms, representing jsp, the 
spin-polarized current density, underline a crucial differ¬ 
ence between the momentum operator and the current 
density operator. These terms can be obtained in several 
ways. They can be derived from a variation of the Hamil¬ 
tonian with respect to A as already done in [47], or they 
can be derived by applying the Foldy-Wouthuysen trans¬ 
formation to the fully relativistic Dirac current j = ec a 
[49[. Here, we adopted a different approach to derive the 
spin-polarized current, see Appendix C. For our purpose, 
however, it is important to note that the additional jsp 
term does not give any contribution. In fact, its matrix 
elements which are required for an ab initio calculation, 
are 

cfi f 

(d'nbspld'm) = ^ J^drV X [4';(r)cr4'™(r)] 

= ^ f dnnx['i>*^{r)eT'P^{r)], (16) 
J s 

with n the unit vector normal to the surface. This in¬ 
tegral vanishes when we consider the infinite surface of 
integration that encloses our sample. Consequently, we 
have proved that, for our purpose, we can work with the 
current operator given by the momentum operator lip. 


D. Optical conductivity tensor and MOKE 

Once the required current density operator, and thus 
the interaction Hamiltonian Eq. (11), is known, linear- 
response theory can be applied to obtain the conductiv¬ 
ity tensor (Ta/ 3 (w). The most accurate way to evaluate 
the conductivity tensor would be to invoke the linear- 
response time-dependent DFT which includes the con¬ 
tribution of the exchange kernel to account for non¬ 
equilibrium electron interaction processes in the excited 
state [50, 51[. However, such calculations on the here- 
desired relativistic level have not yet been performed. 
Also, it is our aim to interpret MOKE experiments where 
the employed laser intensity is less then 10^^ W/m^. In 
this regime the number of electrons removed from be¬ 
low the Fermi energy in Ni is about 0.01 electron [52[ 
and therefore we do not expect that such rearrangement 
of electrons will give rise to a significant contribution 
from excited-state electron-electron interaction. Further¬ 
more, it was shown previously that the bare Kohn-Sham 


linear-response theory described very well the measured 
magneto-optical spectra of metals [46]. 

As it is shown in Appendix B the conductivity ten¬ 
sor within the Kohn-Sham linear-response theory can be 
expressed by an identical equation valid for the fully rel¬ 
ativistic, semirelativistic and nonrelativistic {j = ^p) 
cases. In terms of the Kohn-Sham single-particle energy 
dispersion relations en{k) and matrix elements of the cur¬ 
rent operator jnmi^)y which can be easily obtained in the 
framework of band structure calculations, it reads 

nn'k 

w - lUnn'ik) + i/rl ’ 

where hunn'{k) = Cnik) — e„'(fc). Here the parameter r 
accounts for the lifetime broadening. 

It has already been proven that magneto-optical Kerr 
spectra are well described in a DFT band-structure 
framework using this single-particle formulation of the 
linear-response theory [43, 53[. The polar MOKE spec¬ 
tra are related to the conductivity tensor elements by the 
equation 

$k(w) = 6»k(w) -I- feK(w) 



Here the exchange field is chosen along the local z axis 
and $K is the complex polar Kerr angle that can be di¬ 
vided in the real Kerr rotation 0 k and imaginary Kerr 
ellipticity Ek- Expression (18) is not exact, but in our 
actual calculations below we used the longer exact expres¬ 
sion [46[. Note that in pump-probe magneto-optical ex¬ 
periments the transient MOKE signal is measured, which 
is often taken to be a direct measure of the changing 
atomic magnetization [1, 24[. 

In the following we compute the influence of the rel¬ 
ativistic spin-photon couplings terms on the magneto¬ 
optical spectra of Ni. To evaluate their influence, all the 
calculations described in the next section are performed 
with switching on and off the terms in square brackets of 
Eq. (14). To be precise our ab initio implementation of 
the relativistic momentum operator matrix elements use 
the equation (10) in Ref. [43[ which is exact to all order 
in 1/c^. 

III. MAGNETO-OPTICAL KERR EFFECT 
CALCULATIONS FOR NICKEL 

To elucidate the influence of the relativistic spin- 
photon interaction on the magneto-optical response, we 
preformed ab initio calculations for nickel. In particular 
we compare computed Kerr effect spectra obtained using 
either the nonrelativistic expression of the momentum 
operator p = —ih^ in the evaluation of the conductivity 
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Figure 1. (Color online) Calculated optical conductivity spec¬ 
tra Re[cra: 2 ;(tj)] (top panel) and a;Im[cra;j,(a;)] (bottom panel) 
of fee Ni, using the relativistic and nonrelativistic current 
density formulations. The blue lines are calculated consider¬ 
ing the momentum operator in nonrelativistic approximation 
(p = —iftV). The orange dotted-dashed lines are the calcu¬ 
lations performed using Eq. (14) as momentum operator. 


tensor Eq. (17), which determines the polar Kerr effect, 
or the semirelativistic expression for the momentum op¬ 
erator in Eq. (14). 

It is already well known that the magneto-optical Kerr 
effect is by itself a relativistic effect, since it is directly 
related to spin-orbit coupling [54, 55] present in the ab 
initio calculated electronic structure (particularly, in the 
wavefunctions). The latter we compute with a relativistic 
(four-component) extension of the augmented-spherical 
wave (ASW) code [56], adopting the local spin density 
approximation (LSDA) to the DFT. As has been shown 
previously, using only the nonrelativistic momentum op¬ 
erator in the conductivity tensor calculations (in con¬ 
junction with fully relativistic electronic structure cal¬ 
culations) provides a good description of the MOKE of 
metallic ferromagnets ]46], including Ni. 

In Fig. 1 we show the comparison between the 
interband-only optical conductivity elements, Re[tTa;x(w)] 
and a;Im[(Ta;y (w)], computed with the nonrelativistic mo¬ 
mentum operator as well as with including the relativistic 
corrections to the momentum operator. The calculations 
are performed with a broadening ft/r = 0.03 Ry. As it is 
apparent from the plot the contribution due to the rel¬ 
ativistic terms does not lead to an appreciable change 
in the conductivity spectra. The influence of the addi¬ 
tional light-spin interaction terms on the MOKE spectra 
is shown in Fig. 2. As expected from the results shown 
in Fig. 1, the comparison in Fig. 2 shows that, also for 
the Kerr spectra, the contribution from the relativistic 



Figure 2. (Color online) Calculated magneto-optical Kerr 
rotation and Kerr ellipticity of Ni. The black and blue 
lines are, respectively, the Kerr rotation and ellipticity 
£k calculated without taking into account the additional 
relativistic terms of the momentum operator. The red and 
yellow dashed lines are the Kerr rotation and ellipticity, 
respectively, calculated retaining the relativistic terms of the 
momentum operator. 


correction terms is very small. 

To quantify the influence of the relativistic terms on 
the conductivity tensor we show in Fig. 3 (top panel) 
the difference between the off-diagonal components of 
the conductivity tensor calculated with and without the 
relativistic terms in the momentum operator. This dif¬ 
ference is of the order of 10“^ (xlO^® s“^) for both the 
real and imaginary part. Note that the difference curves 
in Fig. 3 are very smooth and do not show any jitter. 
The reason for its absence is the high numerical accu¬ 
racy in the calculation of Uxyiuj) which is still apprecia¬ 
bly higher (better than 10“^) (for details of the numer¬ 
ical implementation, see ]53]). Hence, there is no doubt 
that we can adequately capture the influence of the rel¬ 
ativistic corrections terms. In Fig. 3 (bottom panel) we 
plot the absolute values of the differences in normal¬ 
ized to the nonrelativistic off-diagonal conductivity, i.e., 
\lS.axy\/\(Jxy nr|) and the same for the complex Kerr angle, 
i.e., |Aci>K|/i$K.nr| = (+ A^)(0| J V A 

For both quantities the normalized differences are of the 
order of 0.1%. Thus, we can conclude that the relativistic 
spin-photon terms do contribute to the magneto-optical 
signal of Ni, but that this contribution is rather small. 
The spin-photon induced change in the MOKE signal 
would consequently be present during the pump pulse, 
where it would be interpreted as as magnetization change 
of the order of 0.1%. 
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Figure 3. (Color online) Top panel: differences between the 
off-diagonal conductivity axyiui) components computed with 
the relativistic and the nonrelativistic current operators. Bot¬ 
tom panel: calculated absolute value of the difference in a^y 
($k) normalized to the quantities \axy,-a-i:\ (|$K,nr|) computed 
with the nonrelativistic momentum operator. 


IV. DISCUSSION AND CONCLUSIONS 

Previous estimations of the influence of the relativis¬ 
tic spin-photon interaction have been made by Vonesch 
and Bigot [39], who considered optical transitions on 
a hydrogen atom within the framework of an extended 
Pauli Hamiltonian. Calculating the matrix elements of 
the spin-photon terms in their Hamiltonian, they found 
that these were of the order of 1 x 10“® (whereas the 
nonrelativistic term p ■ A was of order 1). The largest 
contribution in their treatment originated from the cross 
term Al ■ Aext of the vector potential of the laser radia¬ 
tion and an external magnetic field, which, as mentioned 
before, does not arise in our treatment. Vonesch and 
Bigot estimated a change in the normalized Kerr rota¬ 
tion of 0.5 X 10“^. Thus, even with the additional term 
this estimated change in the Kerr rotation is in overall 
accord with ours for metallic Ni. 

Nonetheless, in spite of the nonzero influence, our ab 
initio calculations do not evidence that relativistic light- 
induced spin-flip transitions could provide a notable de¬ 
magnetization channel. They would appear as a small 
demagnetization effect during the pump pulse which is 
in experiments typically about 70 fs wide. However, dur¬ 
ing and immediately after the pump pulse there will also 
be the influence of “bleaching”, that is, the reduction in 
the optical excitation channels caused by the presence of 
pump-laser excited electrons [52, 57]. The influence of 
such nonequilibrium electron populations on the MOKE 
spectra of Ni have been evaluated previously [58], yet 
without the here-investigated relativistic spin-photon ef¬ 


fects, and were found to be significant. We can thus con¬ 
clude that the nonequilibrium populations have a larger 
effect on the apparent MOKE signal than the relativistic 
spin-photon interaction. 

The demagnetization of Ni after an intensive laser 
pulse has recently been computed by Krieger et al. [59], 
who employed the time-dependent DET formalism. As¬ 
suming extremely intense electromagnetic fields with a 
laser intensity of 10^'^ — 10^® W/m^ they computed an 
appreciably larger demagnetization (of ^ 50 %) than we 
do, which they attribute to dominance of nonlinear ef¬ 
fects. Conversely, in the current investigation we are in 
the moderate fluency regime, with typical laser intensi¬ 
ties of ~ 10^^ W/m^, where the linear interaction Hamil¬ 
tonian Eq. (11) should be sufficient and there is only a 
small number of electrons present in the excited state; 
for this regime our results should be valid. 

Summarizing, we performed the Eoldy-Wouthuysen 
transformation on the Dirac-Kohn-Sham equation in the 
presence of the exchange magnetic field as it is required 
for the relativistic density functional theory in the frame¬ 
work of the local spin density approximation. We ob¬ 
tained a Hamiltonian where several terms are consis¬ 
tent with results derived previously in Ref. [45]. We 
further showed that the spin-polarized term in the cur¬ 
rent density operator is irrelevant for the calculation of 
the conductivity spectra. We discussed the modifica¬ 
tion caused by the relativistic spin-photon terms to the 
linear-response theory for the conductivity, and showed 
that an identical linear-response expression can be ob¬ 
tained for the nonrelativistic, semirelativistic and fully 
relativistic interaction Hamiltonians. We then calculated 
the influence of the relativistic correction terms to the 
magneto-optical Kerr spectra of nickel. In the moderate 
fluency regime, where the linear-response theory is ex¬ 
pected to be valid we find that relativistic spin-photon 
interactions can give a small modification (< 0.1 %) of 
the off-diagonal optical conductivity and of the MOKE 
signal. Thus, our calculations confirm that relativistic 
spin-photon interactions do exist, as originally proposed 
in Ref. 25, but we do not find that these could provide a 
notable channel of laser-induced magnetization loss. 
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Appendix A; The Pauli Hamiltonian 

In the presence of an external electromagnetic field 
(characterized by the vector potential A{r,t)), the fully 













relativistic Dirac Hamiltonian (first without exchange 
field) has the form 

H = ca ■ {p — eA) + (/3 — 1) mc^ + V. (Al) 


In the nonrelativistic limit (which is obtained by perform¬ 
ing the FW transformation), it exactly gives the Pauli 
Hamiltonian 


Tdp = 




consider the semirelativistic case corresponding to the 
extended Pauli Hamiltonian. We also note that, using 
a proper gauge, A = , the above mentioned Hamil¬ 

tonian can be rewritten as first-order interaction Hamil¬ 
tonian in E{r,t). Using the gauge, it is obvious that 
r ■ A = 0, which gives 

-^(r ■A)=0^r-A = r- E 
dV ' 

^ —U ■A = er-E. (B2) 
m 


where the external magnetic held is given by B = V x A. 
If we choose for simplicity a gauge such that 


which fulhlls the Coulomb gauge (V • A = 0) for the uni¬ 
form magnetic held. Then, with this Pauli Hamiltonian 
one can show how the different magnetic contributions 
arise. Namely, the Pauli Hamiltonian can be rewritten 
as 


Therefore, the hrst-order interaction Hamiltonian can 
equally well be written in the form: 

Hi = -e^r,-E = Be-^‘^*\ (B3) 

i 

where are the positions of the electrons and w’*' = u! + 
ijr. In linear-response theory the total average, induced 
current J = jV, with V the volume of the system, is 
computed from (see, e.g. [62]): 


( 2 \ 2 

(A3) 

with g the Lande g-factor, which is 2 for spin degrees 
of freedom. The hrst term obviously is the unperturbed 
Hamiltonian, the dominant perturbation is the paramag¬ 
netic contribution and the last term is the diamagnetic 
contribution [60]. Note that the external magnetic field 
couples to both the spin and orbital angular momentum 
operators, as it should be. 

For magnetic materials the Pauli exclusion principle 
gives in addition rise to magnetic exchange. To include 
magnetic exchange (which dominates over the dipole- 
dipole interaction) the DKS Hamiltonian has to be writ¬ 
ten as [61] 

H = ca-(p - eA) + {l3 - 1) mc^+V+fiBl3 (A4) 

where the exchange field B^^ has to be separated from 
the external magnetic vector potential as it would other¬ 
wise couple to the orbital degrees of freedom. The FW 
transformation of this Hamiltonian leads to the Hamil¬ 
tonian given in Eq. (9). 


Appendix B: Derivation of the optical conductivity 

Introducing the electromagnetic held A(r, t) produced 
by the intensive laser pulse the hrst-order interaction 
Hamiltonian could be written in terms of the momentum 
operator [Eq. (14)] of the unperturbed Hamiltonian, 


Jit) = Tr(po J) + ^ £ dt'{ [Jit), Hjit')] )„, (B4) 


where (...)o means that the average has to be computed 
with the equilibrium density matrix pQ. The hrst term 
refers to the equilibrium current density, which is usually 
taken to be zero in linear-response theory. Note however 
the difference to the derivation given in Refs. [47, 62], 
where this is not done and a second-order interaction 
term A^ is introduced in the Hamiltonian [47]. This 
term is rewritten in Refs. [47, 62] and leads to the Drude 
response (hrst term in (B6) below). Such term should 
however not be included in a linear-response treatment, 
and it is actually not needed, as our derivation shows. 
Our formalism is valid in nonrelativistic, semirelativistic 
and fully relativistic case. We introduce the linear inter¬ 
action Hamiltonian according to (B3) in the second term 
of Eq. (B4) and calculate the integral. Partial integration 
of this second term leads to: 


Jc.it) = -{[J^it),Bit)]) 


0 -iuj+ 


- (B5) 


where B is the derivative of B, which is related to the 
current. Bit) = —e-'^^Tiit) ■ E = —Jit) ■ E. Using Eq. 
(B3) we calculate the commutators and it is evident that 
the integral leads to the current-current correlation in the 
average current, 


Hi = -—U A. 

m 


(Bl) 


Jait) = 


iNe^ 


Eait)e 


-l:fjt'{[J4t),Jpit')]), 


Epit')t 


IbJ^ 


-. (B6) 


This form of the interaction enters in the nonrelativis¬ 
tic, semirelativistic and fully relativistic cases. Here we 
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The conductivity response to the electromagnetic field is Appendix C: Derivation of the spin-polarized 

given as current 


ja{t) = f dt'aap{t-t')Ep{t'). (B7) 

J —OO 


Now, comparing both equations we obtain the linear- 
response expression for the conductivity. Computed in 
Fourier space, the conductivity in terms of (noninteract¬ 
ing) single-particle states is then (see Ref. [46] for details) 


0-a/3(w) 


ie^NSap 

Vmuj~^ 

^ \ " /(Cra) ~ .g ./3 

VhiO+ ^ UJ+-UJn'n 

n.n.' 


(B8) 


where jnn' are the matrix elements of the current den¬ 
sity operator for the single-particle states n and n', 
f{en) is the Fermi-Dirac distribution function of the 
n-th state having energy e„ and huin'n = Cn' ~ fin- 
This linear-response expression is exact for the non- 
relativistic, semirelativistic and fully relativistic cases. In 
the semirelativistic limit, the two terms in Eq. (B8) can 
be approximately joined together, which yields 




hv ^ uj 

nn' ^ ^nn' ^I 

nn 


(B9) 

where the « sign relates to the intraband term i.e., 
n = n') for which the approximation lip Ri p has been 
made. 


The current density operator in semirelativistic form 
was previously shown [48, 49[ to contain a term jsp = 

— ^Cp ^ cr + cr X p). This term also appears in our 
treatment. When we apply the FW transformation to the 
DKS Hamiltonian in Eq. (2), wefindaterm ■^cr-{pxp), 
in the Hamiltonian. This term is taken to be zero for 
obvious reasons. However, it is this term that leads to 
the spin-polarized current density jsp. 

Defining the charge density as p = e5 (r — r) and using 
the Heisenberg equation of motion for the above-given 
Hamiltonian term, 

2 

e.5 (r-f), —a- ■ {p x p) 

g 

= ■ {[S {r - r) ,p\ X p + p X [S {r - r),p\} 

= V ■ {cr X p + ^ xcrj 
2m 

= -V-jsp. (Cl) 

In this derivation, we make use of the fact that, 
the commutator in the position basis, [S{r — r),p] = 
ih\r)'Vr{i'\. Using the continuity equation we extract 
the spin-polarized current density operator as jsp = 

— X p + ^ X cr}. The matrix elements are given 
by" 

/ dri'Sla- X X cr} 

2m Jv 

f>h f 

= ^ J^drV X [4/>4-^]. (C2) 


dp 1 
dt ih 
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